Quantum-critical relativistic magnetotransport in graphene 
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We study the thermal and electric transport of a fluid of interacting Dirac fermions using a 
Boltzmann approach. We include Coulomb interactions, a dilute density of charged impurities and 
the presence of a magnetic field to describe both the static and the low frequency response as a 
function of temperature T and chemical potential /i. In the quantum-critical regime < T we find 
pronounced deviations from Fermi liquid behavior, such as a collective cyclotron resonance with 
an intrinsic, collision-broadened width, and significant enhancements of the Mott and Wiedemann- 
Franz ratio. Some of these results have been anticipated by a relativistic hydrodynamic theory, whose 
precise range of validity and failure at large fields and frequencies we determine. The Boltzmann 
approach allows us to go beyond the hydrodynamic regime, and to quantitatively describe the 
deviations from magnetohydrodynamics, the crossover to disorder dominated Fermi liquid behavior 
at large doping and low temperatures, as well as the crossover to the ballistic regime at high fields. 
Finally, we obtain the full frequency and doping dependence of the single universal conductivity (jq 
which parametrizes the hydrodynamic response. 

PACS numbers; 73.63.-b,05.10.Cc,71.10.-w,81.05.Uw 
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I. INTRODUCTION 

Graphene attracts a lot of interest due to the massless 
Dirac fermions which constitute the low energy quasipar- 
ticles of the undoped materialiii^i^i^. At finite tempera- 
ture and moderate doping, they form a Coulomb interact- 
ing electron-hole plasma of Dirac fermions whose trans- 
port properties are rather peculiar and significantly differ 
from a standard Fermi liquid. Indeed, it has been argued 
that graphene is a quantum critical system^i^i^i^ii^, in 
the sense that the inelastic scattering rate is solely set 
by the temperature and proportional to ksT /h. To date 
the quantum critical and relativistic nature of such Dirac 
liquids, and their dynamic properties, have not been very 
systematically explored, despite the rather fascinating 
possibility it offers to study aspects of relativistic plasma 
physics in an easily accessible tabletop solid state system. 

In a previous work^i, a hydrodynamic approach de- 
scribing the thermoelectric transport properties of such 
systems has been put forward. It predicts a highly con- 
strained form of the frequency-dependent response func- 
tions, and their dependence on doping, magnetic field 
and temperature. The hydrodynamics leaves undeter- 
mined only one transport coefficient in the form of a 
conductivity ctq, whose universal value for an undoped 
system has recently been calculated in Ref. The mag- 
netohydrodynamic analysis further suggests the presence 
of a collective cyclotron resonance whose origin in the 
quantum-critical window is a genuine many-body effect. 
It arises due to particles and holes colliding at a high rate 
with each other while collectively executing an orbiting 
motion at a significantly lower cyclotron frequency. 

The prediction of this interesting relativistic behavior 
was originally initiated by studies of a related relativistic 
systemic which arises in the context of the superfluid- 
insulator transitioUfi^ The putative proximity of such 
a phase transition in the parameter space of various 



strongly correlated systems might significantly influence, 
if not dominate, their low energy physics. This mo- 
tivates the study of linearly dispersing bosonic quasi- 
particles which become massless at the quantum criti- 
cal point. The application of Boltzmann transport the- 
ory to such critical systems was described by Damle and 
Sachde v^^i^^ , and has recently been extended to include 
the effects of a magnetic field by Bhaseen et al.— . With 
small modifications, the analysis presented in this paper 
in the context of graphene can be applied to such sys- 
tems, too. 

The hydrodynamic analysis in Refs . ITllfl^ relied on sev- 
eral assumptions, such as a weak magnetic field, light 
doping and weak disorder, the limits of which remained 
unclear. Further, the conductivity coefficient uq gov- 
erning the entire frequency dependent response was only 
known at the quantum critical point itself, while it was 
not clear how the doping-driven crossover to Fermi liq- 
uid behavior-i2ii^ could be described. These gaps will be 
closed by the present analysis which, similarly to related 
work on graphene and other critical system o-'^^'^'^i^^'^^ , 
starts from a microscopic approach based on the Boltz- 
mann equation, yielding an intuitive physical picture of 
the crossover from quantum-critical to Fermi liquid be- 
havior. 

Not only does the present approach allow us to com- 
pute the coefficient aq and its dependence on chemical 
potential, frequency and magnetic field, but we will also 
determine the precise range of validity of the hydrody- 
namic analysis and the leading corrections to it. The 
Boltzmann approach further allows us to go beyond the 
hydrodynamic regime and to explore the crossover to the 
disorder dominated regime at large doping, where we will 
recover all characteristics of a Fermi liquid. We can also 
study the crossover out of the hydrodynamic regime to 
the regime of strong magnetic fields, which resolves an 
apparent discrepancy between hydrodynamic and Boltz- 
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mann approach concerning the value of the thermal con- 
ductivity of an undoped clean systemi^ii^. 

Several or our results on the transport in the hydro- 
dynamic regime, on the limits of the latter and on the 
behavior beyond hydrodynamics turn out to be very simi- 
lar to exact results which have been obtained for strongly 
coupled, maximally supersymmetric conformal field the- 
ories (CFTs)^^. The latter can be solved exactly thanks 
to the AdS-CFT mapping to a weakly coupled quantum 
gravity problem in a universe that asymptotically be- 
comes an Anti-de Sitter (AdS) space. It is interesting 
that the present weak coupling approach and the exactly 
soluble strongly coupled problem come to very similar 
conclusions regarding the limit of validity of hydrody- 
namics and certain aspects of the physics beyond. 



A. Summary of results 

We briefly highlight our main results that we hope 
to be experimentally observable in the near future at 
high enough temperatures and moderate doping where 
quantum interference effects and disorder are less im- 
portant than the electron-electron interactions. In the 
regime where T > /x, graphene behaves essentially in 
a quantum critical manner and exhibits significant non- 
Fermi liquid behavior, which smoothly crosses over to 
the conventional Fermi liquid physics at larger doping. 
This is reflected in the d.c. conductivity (|6.3p which in 
pure enough samples exhibits a universal, interaction- 
limited conductivity at low doping, and crosses over to 
a conductivity which grows linearly with density upon 
doping. This crossover is predicted to show up in the 
thermopower axx^ too, with significant deviations from 
Mott's law in the quantum critical region, as described by 
Eq. (|6.5p . A similar statement holds for the Wiedemann 
Franz law and the Lorentz ratio. We emphasize that for 
these effects to be seen, samples with rather low disorder 
are needed. Some recent experiments claim that disor- 
der levels for graphene on substrates are currently still 
such that the Dirac physics is dominated by impuritiesi^ . 
However, the recent experimental progress on suspended 
grapheneSS seems a promising route to significantly re- 
duce disorder levels and approaching the regimes where 
the above non-Fermi liquid physics could be observed. 
The latter is expected for clean and large enough samples 
where the inelastic scattering dominates over impurity 
scattering, while the inelastic mean free path (estimated 
in (|5.14p . is still smaller than the sample dimensions. 

Another important prediction of our paper concerns 
the existence of a collective cyclotron resonance in the 
quantum critical regime (in all response functions), which 
smoothly crosses over to the standard semiclassical res- 
onance at high doping. Recent experiments2ii2^i2^ have 
observed a "non-hydrodynamic" cyclotron resonance in 
a regime of strong magnetic fields in which the Landau 
levels and their non-integer spacing in agreement with 
the Dirac equation can be resolved. Our prediction ad- 



dresses however a very different regime at high tempera- 
tures and moderate doping where the quantization of or- 
bits can be neglected. In this collision-dominated, semi- 
classical regime, the collective response of the electron- 
hole plasma averages over the cyclotron frequencies of 
non-interacting particles and holes at typical thermal en- 
ergies. This translates into a resonance frequency pro- 
portional to the doped carrier density, (|7.32p . occurring 
along with an intrinsic, interaction-mediated broaden- 
ing which scales with the square of the magnetic field, 
cf Fig. [71 At high doping the resonance is predicted 
to turn into a sharp peak at the semiclassical value, cf. 
Fig.ini 

Another consequence of the magnetic field is that it 
renders the thermal conductivity k^x finite, even in pure 
systems. At low doping there is an interesting rela- 
tionship between Kxx and the interaction-limited con- 
ductivity, which states that the thermal conductivity 
is inversely proportional to the zero-field conductivity, 
Eq. (|7.I8p . Many of these results can be understood 
from a hydrodynamic point of view. However, we also 
predict crossovers to ballistic regimes when either the 
frequency or the magnetic field are increased such that 
the associated dynamic timescales become shorter than 
the inelastic scattering rate. 

Our paper is structured as follows: In Sec. [TT] we in- 
troduce the model of a Coulomb interacting Dirac liq- 
uid, having graphene in mind in particular. In Sec. IIIII 
we briefly review the derivation of the thermo-electric 
response functions via a magnetohydrodynamic analy- 
sis, with emphasizing the underlying assumptions. The 
formalism of the Boltzmann equation is introduced in 
Sec. IIVI We discuss the two most relevant modes of the 
system, associated with charge, momentum and energy 
currents, and show that due to a peculiarity of 2D sys- 
tems a description restricted to these modes may give 
asymptotically exact results. Sec. |V] discusses transport 
in the absence of magnetic fields and the dependence of 
the inelastic scattering rate on doping. The hydrody- 
namic predictions are recovered and the leading disorder 
corrections are determined. The crossover to the disor- 
der dominated Fermi liquid regime is analyzed in Sec. lVIl 
Finally we discuss the magnetotransport in Sec. IVIIl an- 
alyzing in detail the collective cyclotron resonance and 
the deviations from hydrodynamics at large fields. We 
conclude with a brief summary of the main results. 



II. GRAPHENE WITH LONG-RANGED 
COULOMB INTERACTIONS AND COULOMB 
IMPURITIES 

The effective low-energy description of an undoped 2 
dimensional sheet of graphene is well-known to be cap- 
tured by the Dirac Hamiltonian for massless electrons, 
where the Dirac spinor refers to the pseudospin degrees 
of freedom associated with the two sublattices of the car- 
bon honeycomb lattice. In graphene, the Brillouin zone 
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contains two inequivalent Dirac points ("valleys"), each 
with two spin degrees of freedom, resulting in iV = 4 
species of Dirac fermions. In this paper we consider the 
slightly more general situation of a liquid of weakly inter- 
acting Dirac fermions at a finite chemical potential (dop- 
ing). We assume the electrons and holes to interact via 
standard l/r Coulomb potentials, and allow for the pres- 
ence of charged impurities providing long-range disorder. 
However, we neglect electron phonon scattering which are 
subdominant at low enough temperatures. Their effects 
are discussed, e.g., in Ref. 24. The full Hamiltonian is 
then composed of three parts 



dis ; 



where 



Ho 
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E 



dx 
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(2.1) 



, (2.2) 



with the Fermi velocity vp- The latter was measured to 



be approximately 



21.25.26 



1.1 X 10** cm/s « c/300. 



The spinor representation of the wave-function has the 
following Fourier decomposition 



«'a(x,t) 



Cla(k, <) 
C2a(k,t) 



(2.3) 



where the operators Cia are the electron annihilation op- 
erators on the two different sublattices denoted i = 1,2, 
and a is a multi- index labeling the N fermion species, i.e., 
spin degrees of freedom and the different valleys associ- 
ated with the Dirac points. The formulation of transport 
is simplest in a basis which diagonalizes the Hamiltonian 
Hq. This is accomplished by a unitary transformation 
from the Fourier mode operators {cia, C2a) to the basis of 
chiral particles (7_|_a,7_a): 



cia{k) - -^(7+,(k)+7_,(k)), 



C2a{k) = 



K 
V2k 



(7+a(k)-7-a(k)). 



(2.4) 



We have introduced the following notation: as k is a 
two-dimensional momentum, we can define the complex 
number K by 



K 



ik,. 



where 



{kx ; ky 



(2.5) 



and k — |k| — \K\. Expressing the Hamiltonian Hq in 
terms of 7±aj we obtain 



N 



A=± a=l 



(2^ 



A«FA;7L(k)7,Jk). (2.6) 



In this basis the 1 /r interactions take the for: 

d'^ki d^k2 d^q 
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H, = 



N 

E 

a, 6—1 Ai A2 A3A4 



E 



(27r)2 (27r)2 (27r)2 
xrAiA2A3A4(ki,k2,q)7j;^b(ki -f q)7j;^„(k2 

X7A2a(k2)7Aib(kl) . 



(2.7) 



q) 



Here 



T: 



Ai A2 A3 A4 



(ki,k2,q) 



V'(q,Wki,q) 



1 + A1A4 



1 + A2A3 



X (2.8) 

{K*-Q*)K2- 



|ki +q|/ci 

with Wki,q = w_F(A4|ki -t- q| - Ai|ki|), and 

27re2 



q|fc2 



ere(g,w)|q| 



(2.9) 



is the dynamically screened Coulomb interaction. Here 
er is the dielectric constant due to the adjacent media, 
and e{q,u}) is the dynamic screening function which we 
will discuss below. Note that we have neglected the scat- 
tering between valleys a ^ c, since they involve large 
momentum transfers which are strongly suppressed. 
Finally, we introduce the disorder potential 



dxFdis(x)*t(x)vE-Jx) , (2.10) 



with 



Kiis(x)=^ 



£ X • 



(2.11) 



Here denotes the random positions of charged impuri- 
ties, assumed to be close to the graphene sheet, having a 
charge Ze and average spatial density pimp- Let us also 
express the disorder Hamiltonian -ffdis in terms of the 
7Aa: 



N 



= EEE 



d^ki d^k2 
(2^)2 (2^)' 



■C/AiA2(ki,k2)(2.12) 



i a—1 A1A2 

X exp[ix, • (ki - k2)]7li„(ki)7A2a(k2), 



where 
C^AiA2 (ki,k2) 



2'KZe- 



er|ki 



1 

ka] 2 



1 + A,A2^*^^ 



fclfc2 



(?.13) 



which corresponds to unscreened Coulomb scatterers. 
Note that even though we compute specific results for 
Coulomb interacting particles and Coulomb impurities, 
the formalism easily generalizes to arbitrary isotropic two 
body interactions and disorder potentials coupling to the 
local charge density. 



A. The role of screening 

It is known that in d < 2 generic interactions lead to 
a singularity in the amplitude for coUinear forward scat- 



tering processes 



10.14.27 



To regularize this singularity, we 



will account for screening of the interactions within the 
random phase approximation (RPA) 2^i29 , The dielectric 
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function was calculated in Ref. '29 and has the general 
form 



e{q, Lj) = 1 + a 



'\/ ^Vpq' 



1 f \ ^ ^ 



with /i.2 tending to constants as g 
fine structure constant of graphene: 



0, and a is the 



erhvf 



(2.14) 



which we assume to be small in the present paper, either 
due to a large dielectric constant, or due to its logarith- 
mically small renormalization at low temperature. In the 
static limit and at low temperatures the dielectric func- 
tion reduces to 



e(g,^-0) = l + ^, 

q 



(2.15) 



where qtf — const, x akp ~ an^/^ is the Thomas-Fermi 
wave vector. This leads to screening of charged impu- 
rities at finite doping, replacing the Fourier transform 

of the 1/r interactions by — 27re ^ However, 

since gxF ^ akp, and impurity scattering is dominated 
by momentum transfers of order max[kp, T/vp] the effect 
of screening is small for small a, and will be neglected 
in the following. At larger a, one would however have to 
deal with the screened potential of the charged impuri- 

However, screening plays a crucial role in inelastic for- 
ward scattering processes: In the limit of coUinear for- 
ward scattering where the momentum transfer tends to 
u) — » qvp, the second term in (|2.14p diverges, reducing 
the scattering amplitude to zero. This regularizes the 
logarithmic divergence in the inelastic scattering cross- 
section. Since screening is controlled by the (small) fine 
structure constant a we will neglect it for quantitative 
evaluations that focus on the quantum critical regime, 
and only retain the screening for the purposes of regu- 
larizing the forward scattering amplitude. Having this in 
mind we use the simple approximate form for the RPA 
screened Coulomb interaction: 



^4c(q,t^) 



27re^ 



1 



er|q| 1 ■ 



(2.16) 



where r] (x a will be taken to be an independent small 
parameter throughout this paper. A typical value for the 
proportionality constant 77/0; is 1/16 which is the exact 
result^* for ^ ~ T — 0. Having in mind small a ~ 0.1 we 
will often quote numerical values calculated for a fixed 
value T] = 0.01. 

We note that screening is, however, important in the 
Fermi liquid regime to recover standard Fermi liquid be- 
havior. A comprehensive treatment of this regime, in- 
cluding dynamic screening effects, can be found in Ref.lsil 



III. HYDRODYNAMICS 

We are interested in thermal and electrical transport 
properties of the Dirac liquid subject to interactions, dis- 
order, as well as a perpendicular magnetic field. Each of 
those ingredients scatter electrons out of their linear bal- 
listic motion, and the relative strength of these scattering 
processes defines various physical regimes. We assume 
the external driving force (an electric field or thermal 
gradient) to be applied with a frequency uj which will al- 
ways taken to be small compared to the largest of these 
scattering rates. In particular, we will not be concerned 
with optically driven interband transitions. 



A. Timescales 

The electron-electron interactions induce a finite in- 
elastic scattering rate, which close to zero doping is of 
the order of 



knT 



2KB 



(3.1) 



and thus essentially set by the temperature. This is 
a hallmark of the quantum criticality of the undoped 
graphene system^ti^. At larger doping, when the chemi- 
cal potential ^ exceeds fesT, the inelastic scattering rate 
tends to the familiar Fermi liquid form T~g^ T'^ / ^, if 
the interactions are screened, as is the case in a Fermi 
liquid. As will be discussed in detail in Section fV Al the 
scattering rate is stronger for unscreened interactions. 
However, we will see that only thermal transport is sen- 
sitive to the inelastic scattering rate in this Fermi liquid 
regime, while the electrical and mixed thermo-electrical 
response is dominated by other processes that are deter- 
mined only by elastic scattering from impurities. 

The elastic scattering rate induced by static charged 
impurities is naturally proportional to the density of im- 
purities, and will be shown to be of the order of 



imp 



1 (ZeVe,)Vi.n, 
h max[A;BT, /i] 



(3.2) 



We note that the inelastic scattering rate decreases with 
temperature, while the elastic scattering rate increases. 
The latter is due to the fact that low energy particles are 
more intensely scattered by Coulomb impurities. 

Finally, the "scattering rate" associated with a mag- 
netic field is the typical cyclotron frequency of a ther- 
mally excited carrier. 



,typ 



eB 



max[fcBr, 



(3.3) 



Note that the cyclotron mass in the more familiar expres- 
sion u)c = eB/m is replaced by its relativistic equivalent 
of a typical energy divided by the square of the relevant 
"speed of fight" vp. 
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B. Hydrodynamic regime 



relation 



If the inelastic scattering rate t^^^ dominates, a hy- 
drodynamic description of the low frequency transport 
should apply. This is the case at low enough doping, at 
high temperatures, and in moderate fields. Indeed we 
will show below that the Boltzmann equation recovers 
precisely the predictions of relativistic magnetohydrody- 
namics, with the additional benefit of gaining insight into 
the limits of such a description. In particular, we will find 
that the single transport coefficient ctq left undetermined 
by the hydrodynamic formalism, as reviewed below, is it- 
self a function of the small parameters ujTcc, lu^/^Tcc, and 

^ec/^imp ■ 



C. Response functions from relativistic 
magnetohydrodynamics 

The thermoelectric response of a relativistic fluid in 
the presence of a magnetic field has been derived in 
Refs. [TllfT^ by a magneto-hydrodynamic analysis for the 
low frequency-long wavelength regime. It has been shown 
there^ that the response at small wavevectors k — is 
insensitive to long-range Coulomb interactions, which is 
confirmed by the present analysis. Below we only briefly 
review the important steps and state the results for the 
response functions for future reference. 

Hydrodynamics exploits the fact that on time scales 
much longer than the inelastic scattering time the only 
remaining dynamic modes are the diffusive currents asso- 
ciated with conserved quantities, i.e., charge, energy and 
momentum. The linearization of the conservation laws 
in small deviations from equilibrium then captures the 
diffusive relaxation from long wavelength perturbations 
back to global equilibrium, and allows one to determine 
the frequency dependent thermoelectric response func- 
tions in the hydrodynamic regime 

This program was carried out for a fluid of relativis- 
tic particles (fermions or bosons) in Refs. [Tllfl^ . and is 
briefly reviewed for completeness in App. [X] The deriva- 
tion of the hydrodynamic equations relies on relativistic 
covariance and a constitutive equation expressing the re- 
lation between electrical and thermal currents. In the 
limit of small thermal gradients electro-magnetic fields, 
the specific form of the latter is determined up to a pos- 
itive coefficient ag with units of an electric conductivity, 
cf., (jA6|) . Note that the validity of the relativistic mag- 
netohydrodynamic results is therefore restricted to small 
magnetic fields B. This will be explicitly confirmed by 
the more general Boltzmann theory developed in the fol- 
lowing sections, which will result in a precise criterion for 
the required smallness of B. 

The thermoelectric transport coefficients describing 
the current (J) and heat current (Q) response to elec- 
tric fields and temperature gradients are defined by the 



a a 
Ta K 



E 



(3.4) 



where a, a and k are 2x2 matrices acting on the spatial 
X, y-components of the driving fields. Rotational invari- 
ance in the plane imposes the form 



1 



(3.5) 



where 1 is the identity, and e is the antisymmetric tensor 
with e^jy = ^^yx ~ 1- Note that the two off-diagonal 



entries in (j3.4p are related due to Onsager reciprocity. 

The thermal conductivity, k, defined as the heat cur- 
rent response to — VT in the absence of an electric current 
(electrically isolated boundaries), is given by 



k = K — Taa ^a. 



(3.6) 



The analysis of the hydrodynamic equations in a weak 
magnetic field yields the following frequency dependent 
response^ii— The longitudinal and Hall conductivity are 



CTxyi^) 



u {uj + 1^ + ii^lh) 
{uj + i-if - ujI 
^ — 2i^u) 



P 



B (^uj + i'-f) — cj^ 



(3.7) 
(3.8) 



The thermopower and transverse Peltier coefficient are 
found to be 

f ^ _ ^ [o"q(m/^)(^ + h) - isp/je + P)] 

OLxxK^J) — ; —2 ,{6.\)) 



axy{^^) = - 



S ujl+-i'^ 



- fip/{sT)] 



and the thermal conductivities 

,£ + P • S^T 



(3.10) 



Kxxi^^) = 

K-xyii^) = 
Kxxii^) = 



-7^ - ij^ ^ + gq^c^ (c^ + Z7) 
{uj + ijf - Lol 



T (tj + i 7) - Wc 

. {e + P) {uj + iiolh) 

je + P) LOc 



(3.11) 
(3.12) 



In these formulae a collective cyclotron frequency ujc and 
a damping rate 7 have been defined as 



eBpvp 

¥+p) 



7 ^ 



{e + P) 



(3.13) 



where here B is given in SI units. To obtain results in 
more customary cgs units one should replace B B/c 
throughout the paper. This shows that the speed of light 
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merely plays the role of a coupling constant determin- 
ing the strength of the magnetic field. By introducing a 
phenomenological relaxation rate t-^^ (due to weak im- 
purity scattering) into the momentum conservation law, 
one finds that the above formulae are simply changed by 
the replacement lu — s- oj-f i/rimp. However, since disorder 
breaks explicitly the relativistic invariance by singling out 
its own rest frame, it is to be expected that the resulting 
expressions for the response are merely qualitatively cor- 
rect and become exact only in the limit when the elastic 
scattering rate is the smallest rate in the problem. This 
will indeed be confirmed below. 

It will be one of our aims in the subsequent sections 
to rederive these response functions in the appropriate 
hydrodynamic regime, and to establish the precise lim- 
its of validity of the hydrodynamic description, and in 
particular the admissible range of B and to. Further we 
will obtain explicit expressions for the transport coeffi- 
cient ctq as a function of the system parameters, most 
importantly as a function of the chemical potential /i. 
This will shed light on the crossover from relativistic, 
quantum critical response in the regime ^ T where 
both particles and holes contribute to transport - to the 
Fermi liquid regime |/^| 3> T where only one kind of par- 
ticles contributes. In the latter regime we will recover 
the standard laws governing Fermi liquids. We will also 
study the disorder dependence of the response functions 
and discuss the crossover from the interaction dominated 
to the disorder dominated regime. 



IV. BOLTZMANN TRANSPORT 

A. Applicability of Boltzmann transport theory 

For the Boltzmann equation to be valid, one requires 
the existence of well-defined, sharp quasiparticle exci- 
tations and sufficiently weak interactions so that scat- 
tering does not lead to strong many body correlations. 
This framework is of course much more general than 
hydrodynamics. We therefore can extend our study of 
transport into regimes of strong magnetic fields and dis- 
order. However, we are always restricted to a regime 
where fc^T is much bigger than the cyclotron energy of 
thermal particles, i.e., ksT ^ huj^^. This ensures that 
we do not need to account for Landau quantization of 
electron orbits and the quantum Hall effect, which lies 
beyond the semiclassical Boltzmann equation. Similarly, 
localization corrections which derive from quantum in- 
terference cannot be captured by a simple Boltzmann 
approach. However, as long as the interference effects do 
not drive the system insulating, as ascertained in several 
recent studies of non-interacting electrons, at least for 
random point-like disorder'^^iSS:, one can expect that the 
Boltzmann approach applied to a field theory with ap- 
propriately renormalized parameters for interactions and 
disorderi^SiSl will capture a large part of the phase dia- 
gram in temperature, disorder and interaction strength. 
Such an approach was recently taken to predict a loga- 



rithmic increase with temperature of the conductivity in 
clean undoped graphene^". 

The central object in Boltzmann transport theory is 
the distribution matrix of the quasiparticles 

/AA'(k,t) = (7L(k,t)7,,,(k,i)) , (4.1) 

where there is no sum over a on the RHS, and we assume 
the distribution functions to be the same for all valleys 
and spins. For all further discussions we will neglect the 
matrix elements off-diagonal with respect to the helicity 
basis labeled by A, which brings us back to the familiar 
quasiparticle distribution function 

A(k,0 = (7ljk,t)7,Jk,t)). (4.2) 

This approximation can safely be made since we are 
interested in low frequencies fUv <^ ksT where field- 
induced, coherent interband transitions leading to off- 
diagonal correlations (with f\^\') can be neglected. In 
equilibrium, i.e., in the absence of external perturba- 
tions, the distribution functions are Fermi functions 



/A(k,i) - fl{k) 



(4.3) 



at the finite chemical potential /i, as defined by the dop- 
ing or gate potential. 

We consider the Boltzmann equation in the presence 
of an electrical field E, the Lorentz force due to a perpen- 
dicular magnetic field B — Be^, and a spatially varying 
temperature T(r): 



^ + e (E + VA,k X B) ~ + VA,k • |: ) /A(r, k, t) 



-Xcon[A,r,k,i|{/}]. 



(4.4) 



Here Xcoii[A, r, k, t \ {f}] denotes the collision integral due 
to Coulomb interactions and impurity scattering, and 
VA,k — Vk^Ak = Xvpe^ with ek — (k/fc), denotes the 
quasiparticle velocity. We rewrite the equation specify- 
ing the driving gradients on the RHS: 



d 

dtf\{r,'k,t) + eXvpB {ei^ x e^) • — /A(r,k,t) 
+Xcou[A, r, k, t\{f}] = • eE + • VT, (4.5) 

where 

= A^ek/,"(fc)[l-/,"(fc)], (4.6) 

and 

= - ^"^%V^ ek/A°(fc)[l-/A°(fc)]. (4.7) 

We seek to solve (|4.4p in linear response and thus param- 
eterize the deviation of /A(r, k, f) from its equilibrium 
value in the standard way"^* as 
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/,(r,k,c^) = 2^5H/«(fc,r(r)) + /°,[l-/0,]^ek 



eE(w)g 



(E) {vpk 
I.A I T 



vr(w) 



/A°.[l-/A\]^(ekxe,). 



(4.8) 



with r(r) = T + r • WT{uj) and dimcnsionlcss functions 
5a(A). 

We note that exactly at particle hole symmetry (/i = 

0), an applied electric field generates perturbations 
having opposite sign for quasiparticles and quasiholes, 



(E) 

9x 



vpk 



(E) ( Hfk 



(4.9) 



whereas a thermal gradient will generate symmetric per- 
turbations, 



(4.10) 



However, in the case of a finite chemical potential the 
distribution function will have a generic dependence on 
A. Notice also that the perpendicular components of the 
perturbations, g±^\, vanish in the absence of a magnetic 
field. 



B. Matrix formalism 

In this section we will set up the calculational frame- 
work for all subsequent discussions. A standard way to 
deal with integro-differential equations consists in ex- 
panding the solution gx into a set of basis functions 

</>n(A, k), 

n 

g±Ak) - ;^^^0„(A,fc) , (4.11) 



and to express the integral equation as a matrix equa- 
tion, by multiplying it from the left with different basis 
functions and integrating and summing over k and A, 
respectively. 

From now on we take momenta k to be given in units 
of ^ unless stated otherwise. More generally, we will 
use units in which h = ksT = vp = ^, but restore those 
in final results. The scattering terms then turn into a 
matrix acting in the space of expansion coefficients ipH 
and V'n J which we organize into a doublet of vectors 



(4.12) 



allowing us to cast the Boltzmann equation into the com- 
pact form: 



M -BY rE,T 
B M j ^ 



yE.T 





(4.13) 



We have used that in linear response we can deal with 
the response to the electric field and the thermal gradient 
separately. A matrix inversion yields the solution 



'E.T 



i - 




{'7) 


\ -K 


K ; 





(4.14) 



with 



K = (M + BM-^B) ^ , 

K = M-^B{M + BM-'^ByY (4.15) 

In the above, the matrix A4 is the sum of three terms: 

M = M^'^ + X™P + M"^, (4.16) 

the first accounting for inelastic scattering due to 
Coulomb interactions, the second describing impurity 
scattering, while the last derives from the time deriva- 
tive in the Boltzmann equation. 

Let us consider the first two terms, corresponding to 
the collision term on the right hand side of (|4.4p . For 
sufficiently weak interactions and for dilute enough im- 
purities, the collision integral is given by an application 
of Fermi's golden rule to two body collisions, and to the 
scattering from Coulomb impuritieai^. The correspond- 
ing matrix elements are given in App. |b1 The strength of 
the inelastic scattering rate due to electron-electron in- 
teractions is characterized by while elastic scattering 
rate from impurities is measured by the dimensionless 
parameter 



(4.17) 



The magnetic field deflects particles from linear prop- 
agation at a rate proportional to the dimensionless pa- 
rameter characterizing the magnetic field strength 



eBvl 



(4.18) 



The relative magnitude of , A and b defines various 
transport regimes which we will discuss below. 
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C. Linear response 

The heat current (Q) is related to the energy current 
(J^) and the electrical current (J) via 



Q = 

e 



(4.19) 



Given a perturbation of the distribution function 
parametrized by the associated heat and electrical 

currents are given by the expressions 



Qii,i- 



^E/|^S/.Vl-A\) (4.20) 



where the shorthand 



fO ._ 
Jxk ■ — 



(4.21) 



denotes the Fermi distribution. 

Using the matrix elements of the driving terms, cf., 
(|B10[) . we can express this as 



J|l,-L 



""5^ ), (4.22) 



Using the formal solution (I4.14|) for the coefficients -01 
we immediately read off the longitudinal transport coef- 
ficients defined in 



J.,{u) NekB ^ 



Kxxiuj) 



VxT{uj) 2h 
Qx{uj) 



NklT - 



-Tt-KTt- (4.23) 



In the presence of a magnetic field, the transverse trans- 
port coefficients are finite as well and given by 

axy(uj) = ^ , , = -^-Fe ■ K--Fe 



Ey{Lj) 2h 



Jx{uj) 


Neki 


VyT{u) 


2h 




Nk%T 


WyT{cu) 


2h 



Tt ■ KTt ■ (4.24) 



D. Choice of basis 



In order to analyze the response functions, it proves 
essential to choose a well-adapted basis 4>n=Q j\f{^,k) 
to expand g\\,±-x{k) into. 



The structure of the currents and of the driving terms 
(I4.22P suggests to use the modes 



bo{X, k) = fc, 
^i(A,fc) = A. 



(4.25) 
(4.26) 



Moreover, it will be convenient to complete the basis in 
such a way that the (j)n>2 do not contribute to the electri- 
cal and thermal currents. Due to reciprocity, this implies 
in turn, that these modes do not couple to the driving 



fields, i.e., ^^^2 — O'" 




(4.27) 



According to (|B10[) this basis choice imposes the fol- 
lowing two constraints for all n > 2 

E / d^fc/Afc(l-/Afc)</'n>2(A,fc)0o,i(A,fc) = 0.(4.28) 

A •' 

From the expressions in App. [B] it is easy to check that 
this implies the vanishing of the matrix elements 



MZn = MZ, = 0, 
Bn.n = 0, for aU n>2. 



(4.29) 



This will play a vital role in establishing the relativistic 
magnetohydrodynamics in Sec. IVIII 

Note that the part of the collision integral due to elastic 
scattering from disorder does not follow the same pattern. 
Rather it has non- vanishing matrix elements A^g™'' for all 
n. This is closely related to the breaking of relativistic 
invariance by the disorder. 

The mode (j)o{X^k) is central not only because it nat- 
urally describes the energy current, but it also has the 
further important property that it is the deviation gen- 
erated by a transformation to a moving frame, and trans- 
lational invariance protects this mode from decaying due 
to Coulomb interactions. Indeed, one immediately checks 
that momentum conservation implies that (/>o is a zero 
mode of the Coulomb collision operator, see Eq. (jB3p . 



(4.30) 



However, the momentum is not conserved by impurity 
scattering and the magnetic field. 

The relativistic linear dispersion is essential to ensure 
that the "momentum mode" 0o coincides with the "en- 
ergy mode" entering the energy and heat currents. More- 
over, it is because of this relativistic dispersion that the 
scattering due to the magnetic field couples (po only to 
the "electrical current mode" This structure in the 
Boltzmann equation is crucial to retrieve the relativistic 
magnetohydrodynamics, as will become clear in Sec. lVIll 



E. Leading logarithmic approximation 

A significant simplification occurs in two dimensions, 
as noticed in a preceding publication^". Indeed, a close 
analysis of the two body collision integral (|B3|) shows that 
there is a logarithmic divergence in the phase space for 
nearly coUinear forward scattering processes^"'^'* which 
is only cut off at small angles of order 0{a), e.g., by 
the RPA screening of the Coulomb interactions, see 
Eq. (|2.16p . Interaction corrections to the linear disper- 
sion or the inclusion of a finite life time of the quasipar- 
ticles would provide a similar cut-off, too. In the limit of 
weak interactions a <C 1, but still in the hydrodynamic 
regime (assuming even smaller impurity scattering rates, 
frequencies and magnetic fields, A,uj/T,b ^ a^), this 
leads to an equilibration among excitations moving in 
the same spatial direction. The off-equilibrium distri- 
bution function will then be characterized by an angle- 
dependent effective temperature and chemical potential. 
Since the angular dependence in linear response has to be 
proportional to the projection of k onto the driving field, 
such perturbations to the distribution functions corre- 
spond precisely to the two modes (po and ipi introduced 
in the previous sectio n Provided the logarithmic 
anomaly in the forward scattering is sufhciently strong 
the above justifies to restrict the analysis of the Boltz- 
mann equation to these two modes. 

To be more precise, we actually invoke that the in- 
elastic relaxation of all other modes (j)n>2 is faster by a 
factor of the order of log(l/77) than that of the "electrical 
current mode" 0i, where cx a is the screening param- 
eter introduced in (|2.16p . Up to logarithmically small 
admixtures of other modes, corresponds indeed to the 
eigenvector of A4'~^^ with the smallest positive eigenvalue. 

It follows that we may obtain exact solutions of the 
Boltzmann equation to leading order in [log(l/?7)]^^ by 
restricting ourselves to the Ansatz 




Figure 1: Dependence on the screening parameter -q of the 
conductivity of undoped graphene, a^a{fj, = 0; r;) in units 
of e'^/h. At very small rj, a a{0; i]) approaches its limiting 
value^" 0.121 as a^[cr(0; 77) - a(0; 0)] - [log(l/ry)]-\ the solid 
line being a fit to a + 6[log(l/r;)]~^. The lower curve is a linear 
fit to the data obtained from the two mode approximation, 
which becomes asymptotically exact as 0. 



boiling down to using the single mode (j>i. We will 
thus use the Ansatz (|4.3ip later to describe analytically 
the crossover from interaction-dominated to disorder- 
dominated regime. However, for most explicit results, 
and for the entire analytical discussion of magnetotrans- 
port we will not resort to this approximation. Never- 
theless, it will become clear that the essential physics is 
captured by the dynamics of the two modes (jjo , (t>i which 
carry the information of energy and charge currents, re- 
spectively. The matrix elements pertaining to these two 
modes are given explicitly in App. [B]and will be used in 
the discussion of transport in the remaining sections. 



= VoMfc + V^iHA, (4.31) 



TRANSPORT IN THE ABSENCE OF A 
FIELD 



for the deviation of the distribution function. This is 
demonstrated in Fig. [1] where we plot a'^a{fi = 0) as a 
function of 77. The lower data set was evaluated within 
the two mode approximation as in Ref. [Tol . while the up- 
per curve was obtained by solving the Boltzmann equa- 
tion projecting onto 12 basis functions (/>,i(A, fc) and in- 
verting the resulting matrices 4^ It is clearly seen that in 
the limit of weak screening 77 ^ 1, the two mode approx- 
imation becomes exact. 

This two mode approximation is well justified in the 
hydrodynamic regime where the inelastic scattering dom- 
inates the dynamics. However, we will argue in Sec- 
tion IVI Al that the Ansatz (|4.3ip also captures the re- 
laxation time approximation which is widely used in the 
disorder dominated Fermi liquid regime /i 3> T. In 
the case where the inelastic scattering can be neglected, 
the relaxation time approximation even becomes exact. 



As we will see in later sections, it is essentially the slow 
dynamics of the momentum mode (j>o which gives rise to 
the special relativistic structure of the response functions 
in the hydrodynamic regime. It is thus convenient to 
treat the zero mode 0o separately, and write the matrix 
M ^ M^^ + M^'^ + 7W'"^P in the form 



M 



( Mo Ml M2 

Ml Mil M12 

M2 M21 M22 

Vm„ M„i M„2 



Mi„ 

M2„ 



(5.1) 



Owing to (|4.27p . we only need to know the inverse of this 
matrix in the sector spanned by (/)o , (j)i in order to calcu- 
late the response functions (|4.23l4.24p . In this sector it 
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assumes the form 



M- 



/ 1 a * 

a c * 



with 



= -^A/„(M-i)i„ 

n>l 

m,n>2 



(5.2) 



det[M(i"'i'")] 



detM 



where det[M(^"'^™^] are the subdeterminants of the ma- 
trix M where the two rows 1 and n, and the two columns 
1 and m have been dropped. The above expressions are 
very useful in the hydrodynamic limit, where typical ma- 
trix elements of M are much larger than the entries Af„ . 



A. Clean case 

We recall that in the absence of disorder all M„>2 van- 
ish, which simplifies the above expression for the matrix 
inverse to 



91 



1 



I 1 -Migi * 



Mo - Mlgx 



(5.3) 



Note that in this case the judicious choice of the ba- 
sis 4)n allows us to summarize the effect of all the modes 
4'n>2 into a single (frequency dependent) matrix element 
(7i, which enters all the response functions. 171 character- 
izes the inelastic scattering rate due to electron-electron 
interactions as we will detail below. In the case where 
the inelastic scattering rate dominates, it is convenient 
to write gi = 51 / to exhibit explicitly the scaling with 
interaction strength, gi being a number of order 0(1). 

In general it is sufficient to use a relatively small num- 
ber of bases to obtain a good accuracy for gi. However, 
as we discussed in Section IIVEI in the hydrodynamic 
regime in two dimensions, there is even a further simpli- 
fication that allows us to neglect the remaining modes 
4>n>2 to leading order in log(l/77). 

From Eqs. (|4?T4l) . (|BI0)) . and (|4?23l) it is easy to show 
that the response in a clean system is given by 



c(a;;/Lt, A = 0) = e 



2 P^«F 1 



(5.4) 



where aq = aQ{^,,uj) is the lo- and /i-dependent trans- 
port coefficient 



o'q(m, 



1 



1 — iuJTp 



-, (5.5) 



which was left undetermined in the hydrodynamic for- 
malism, cf. p.7p . In the above we have defined the in- 
elastic scattering rate as 



_i 2 N ksT 



251 fi 



iVln[2cosh(/i/2)] p'^[hv) 



27r 



{e + P)T 



(5.6) 



B. Inelastic scattering rate in the Fermi liquid 
regime jj. ^ T 

In the above expressions, gi is a scaling function of 
fi/T and wToo. As long as we concentrate on the hydro- 
dynamic frequency regime, the latter is negligible. How- 
ever, the /i-dependence of gi and Toe are important. 

An order of magnitude estimate for gi can be obtained 
from the inverse of the expression (IB12|) for Mf;f'. If the 
electron-electron interactions are not screened (as is the 
case to lowest order in a), this multiple integral satu- 
rates to a finite value at large fJ-/T. This is explained in 
detail in App. [Cl and is borne out by direct numerical 
evaluation, 



gri(M»r) = o(i). 



(5.7) 



The inelastic scattering rate for unscreened long range 
interactions can similarly be estimated to scale as 



(5.8) 



as follows from the above together with (|5.6p . 

If we include screening of the interactions one finds 
instead 



and the analogous estimate for t^^^ 
Fermi liquid behavior 



(5.9) 



'cc,sc 



yields the familiar 



(5.10) 



We have calculated the full function g^^{fj,) in the limit 
of vanishing screening (i.e., 77 ^ 0), where the evalua- 
tion via the formula for the matrix element (|B12[) be- 
comes exact. The resulting static transport coefficient 
ctq(w — 0, /i) is plotted as a function of fi/T in Fig. [51 
It is interesting to note that in the regime fi ^ T, oq 
decays as {T/ ^mY, reflecting that relativistic physics and 
quantum criticality associated with the presence of par- 
ticles and antiparticles becomes less and less relevant. 

The above discussion refers to the typical relaxation 
time for modes g{k) which vary significantly over a range 
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Figure 2: The normalized transport coefBcient aQ{iJ.,u! = 0) 
as a function of n/T. 



of order T around k = fi. These modes are associated 
with so-called "vertical processes" that degrade thermal 
currents very efficiently^. However, their contribution to 
the electrical current is small, which is reflected by the 
smallness of ctq scaling as (T///)^. The electrical con- 
ductivity is thus largely dominated by the ffist term in 
(|5.4p . and this holds true also in the presence of weak 
disorder, see Eq. (|5.16p below. This statement is in- 
dependent of whether electron-electron interactions are 
screened. However, the same is not true for the ther- 
mal conductivity which turns out to be sensitive to un- 
screened Coulomb interactions even when /i 3> T. This 
will be discussed further in Sec. IVIDI 



C. Quantum critical regime 



where we have used gi{0) « 1.24. Accordingly, the in- 
elastic mean free path evaluates to 



2.3 



2r[K] 



/im . 



(5.14) 



D. Weak disorder 

In a clean system, the denominator in the expres- 
sion (|5.3p for the conductivity vanishes as a; ^ since 
Mq,Mi ~ Lo. This reflects momentum conservation due 
to which the mode cjiQ does not decay. 

However, the presence of disorder shifts the pole from 
uj — into the negative halfplane, such that the denomi- 
nator behaves as a; -I- zTj" p. This defines the elastic scat- 
tering rate T~^p. To leading order in weak disorder it is 
given by 



A{p++p-) Pi, 



p h 



0+ +p- 
e + P 



max[/cBT, p] 



(5.15) 



One can explicitly verify that in an expansion in weak 
disorder, i.e., in T^p/T~g^ ^ A/a^, the conductivity is 
given by 



P vj. 



T- * — lUJ £ + P 
' imp 



(Tq +5a{A,uj,p) , 
(5.16) 



Note that close to particle-hole symmetry, p,<^T, t^^^ 
is essentially given by the temperature, r^"^ ^ a^kBT/h, 
cf. Eq. (|5.6p . This is a hallmark of many quantum critical 
systems where often a relativistic effective theory emerges 
(as characterized by a dynamical critical exponent z — 
l)ii with the temperature T being the only energy scale 
left. 

At particle-hole symmetry, the above expressions 
(I5.4l5.5p become identical to the results reported in a 
previous publication^*^ for the clean limit: 

a^^{uj,p.^O) ^ a^- — g-. (5.11) 

one finds the explicit numerical values 



In the limit rj 



1 N , 
0.121 e2 



log(2) 

TT 

_ 0.760 e2 

h h ^ 

h ffi(0)log(2) 



0.274 



(5.12) 
, (5.13) 



with 6cr{A,oj,fi) =0{A/a^). 

Similarly one finds for the other thermo-electric re- 
sponse functions 



axx{uj;p,A) 
Krcx{uj;p,A) 



spvj, 



imp 



,e + P 
1 sWp 



imp 



e 1 



P 



+ 5k{A,uj,p), 



with disorder corrections of order 5a, Sn — 0{A/a'^). 
Remarkably, by dropping these higher order terms, we 
recover precisely the expressions predicted by relativis- 
tic hydrodynamics with a phenomenological momentum 
relaxation rate implemented via uj ~^ lj + ifl^^. Since 
disorder breaks the relativistic invariance of the particle- 
hole plasma, it is not surprising that it eventually spoils 
the relativistic structure of the response functions when 
the disorder-induced elastic scattering rate t~^^ becomes 
comparable to . 
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VI. CROSSOVER FROM QUANTUM CRITICAL 
TO 'FERMI-LIQUID' REGIME AT > T 

A. Two mode approximation 

As mentioned before, in the Boltzmann approach we 
are not restricted to smah disorder satisfying A/a^ ^ 1. 
In order to obtain insight into the crossover from the 
hydrodynamic to the disorder dominated regime, we will 
adopt here the two mode approximation. The latter is 
not only convenient for analytical treatments, but also 
turns out to describe both the interaction and disorder 
dominated limit well. 

In a regime dominated by elastic scattering one usually 
resorts to the relaxation time approximation. 



gx{k) ^ -{df/dE)TkVkeE. 



(6.1) 



In the case of unscreened Coulomb scatterers this is exact 
if the relaxation time is chosen to increase linearly with 
the energy, Tk ~ k. In the 'non-relativistic' limit jj, 3> 
T where thermal excitations of only electrons or holes 
are relevant, this is essentially equivalent to making an 
Ansatz with the single mode (jji, and the corresponding 
relaxation time is related to the coefficient tpi in (|4.1ip 
by 



1^] 



I hvp k H 
knT kuT 




(6.2) 



10 Pimp 



Figure 3: The d.c. conductivity, o?axx in units of e^/?i, as a 
function of the ratio of doped carrier density to impurity den- 
sity. The disorder strength and temperature have been chosen 
such that A/of^ = 0.25. The plot shows the crossover from 
the quantum critical regime, where the conductivity is essen- 
tially limited by inelastic scattering between electrons, to the 
regime dominated by elastic scattering from Coulomb impu- 
rities at higher doping. In the latter regime the conductivity 
increases linearly with doped carrier density, see Eq. (|6.4|) . 
The red and yellow curves correspond to the two terms in 
Eq. (|6.3|l . respectively: the contribution from modes relaxing 
due to inelastic scattering, and the contribution from the mo- 
mentum mode, which is only disorder limited and dominates 
at large density. 



On the other hand, in the regime where electron- 
electron scattering dominates we have already argued 
that enhanced forward scattering allows us to restrict to 
the two modes (/)o,i to leading order in l/log(7y). This 
suggests that an approximate description of the Boltz- 
mann equation (|4.4p restricted to those two modes ac- 
tually yields a rather accurate description of the whole 
crossover from quantum critical to disorder dominated 
regime. 

In this approximation the full expression for the longi- 
tudinal conductivity assumes the form: 



1 



P Vp 



T- ^ — lUJ £ 
imp 



p 



(6.3) 



A/<^'-iw(e+P) 



[A/i^'-iwp] 



the first term diverges in the clean static limit except 
at the particle-hole symmetric point where its numerator 
vanishes. 

At zero doping the above two-mode approximation re- 
duces to the same expression as (|5.1ip . with the replace- 
ment a^/gi ct^/gi + ^/(27r). This means that an 
impurity scattering rate proportional to A is added to 
the inelastic scattering rate in the denominator. The 
purely interaction dominated, 'quantum-critical' con- 
ductivity is thus visible only in weak disorder where 
(A/Q;^)gi/(27r) <^ 1; otherwise disorder dominates the 
response at all dopings. The crossover from interaction 
dominated to disorder dominated transport occurs when 
the two terms in (|6.3p are approximately equal, i.e., when 
P ^ Pimp, cf. Fig. H 



B. Linear conductivity at large doping 



with Tj^p as given in (|5.15p and the functions I^^\fi) 
defined in App. [B) In the second term we have set the 
factors oi h = vp — ksT ~ 1. We have decomposed the 
response function into two parts: The first term is inde- 
pendent of the inelastic scattering rate (that is, of 51), 
and is thus entirely determined by impurity scattering. 
The second term, <t^^, has a finite d.c. limit as A — > 0, 
and it reduces to the quantum critical value of conduc- 
tance of pure undoped graphene as p, — > 0. However, 



We can make direct contact with earlier studies which 
considered disorder dominated transport in the non- 
interacting Fermi liquid regime of graphene. From 
Eq. (|6.3[) and Eq. (jB20p it is easy to check that at large 

^ T, the first term in Eq. ()6.3|) dominates the static 
conductivity which is thus entirely determined by impu- 
rity scattering. The full conductivity and the contribu- 
tions of either term are plotted in Fig. [3l 

In the disorder dominated limit we recover the expres- 
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sion for the conductivity of non-interacting electrons in 
doped graphene in the presence of Coulomb impurities, 
previously reported in Ref. 



2 2 2 

e + P 



TT (Za)^ h Pimp 



(6.4) 



Note that in our approximation where we neglect screen- 
ing of the impurities, there is no constant offset to the 
term linear in p. It would be straightforward to include 
RPA screening in the impurity potential, which is known 
to modify the numerical prefactor of the linear density 
dependencei^ii^i^ in (|6.4p and produce a positive offset 
proportional a^, see Ref4^ 



3 Oxx 




T 



C. Mott's law 

In the regime /i 3> T we expect to recover ordinary 
Fermi liquid behavior. In the approximation with two 
modes ipo,!^ thermopower Uxx is given by 



axx{fJ-,uj = 0) = 



e + P 



(6.5) 




where in the second term we have again dropped fac- 
tors of h,VF,kB and T. At large fi^ T, the first term 
dominates again, similarly as for the conductivity. 

In a Fermi liquid where scattering at low temperature 
is dominated by impurities, one expects the thermopower 
to be related to the conductivity according to Mott's law. 



,{fi,LU = 0) = -—k%T 
3e 



da{ii, w = 0) 



(6.6) 



Indeed, using Sommerfeld expansions for p and s in 
the above results, one can easily show that this rela- 
tion holds for /i ^ T. Not unexpectedly, however it 
fails in the quantum critical regime fj, < T. The ra- 



tio 



is plotted in Fig. 2] as a function 
of fj,/T (assuming unscreened electron-electron interac- 
tions). Note that a^x tends to zero as fj, vanishes. The 
same holds for da/dfi{fi) but the ratio of the two quan- 
tities remains finite. 



Figure 4: The ratio R 



as a function of 



da / d^(fj,) 

fj,/T. In the Fermi liquid regime (^ 3> T) -R tends to 1 as 
predicted by Mott's law. 



diverges as /i ^ instead of acquiring the standard Fermi 

liquid value of ^^fr-^- This is another peculiarity of the 
quantum critical regime. 

One would expect to recover the Fermi liquid value 
of the Lorentz ratio, or in other words, the Wiedemann 
Franz law, by going to sufhciently low temperatures in 
the Fermi liquid regime p 3> T. This is indeed so if we 
assume the electron-electron interactions to be screened, 
such that inelastic scattering processes degrading the en- 
ergy current become inefficient as compared to impurity 
scattering at low T. This is also in qualitative agreement 
with the thorough analysis of Ref. 31, which in addition 
discusses the contribution of bosonic particle-hole exci- 
tations to thermal transport. 

If the interactions are not screened, however, the 
Lorentz ratio remains sensitive to them and does not 
reduce to the Fermi liquid value. The reason is as fol- 
lows. Even though an expression analogous to (|6.3I6.5|) 
exists for 'Kxx, and is dominated by its first term 
s^Tw|,ri,„p/ (e + P), the same is not true for the thermal 
conductivity Kxx — ^xx ~ Ta'^x/'^xx which is sensitive 
to the subleading corrections to the dominant terms in 
(JxxTCtxx and Kxx- If the electron-electron interactions are 
not screened, implying that hT~J^{p) ~ /i for large /x, the 
effects of these interactions persist in these subleading 
terms even deep in the Fermi liquid regime ^ T. 



D. Wiedemann- Franz law 



VII. MAGNETOTRANSPORT 



In the clean limit, the Lorentz ratio 



L : 



(6.7) 



In this section we turn to transport in the presence 
of a magnetic field. The matrix formalism with the ap- 
propriate choice of basis functions will prove particularly 
useful to recover the hydrodynamic response and to find 
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explicit expressions for the more general structure of the 
thermo-electric response functions. To this end, we need 
to develop the matrix formalism of SecHvl We anticipate 
that the dynamics of the momentum mode (j)Q captures 
most of the physics in the hydrodynamic regime. It is 
therefore convenient to introduce a notation that singles 
out the components pertaining to it: 



l|:-L 



(7.1) 



where D ~ E,T. We recall that our basis choice implies 
= {F[',0, . . . , 0)^. Further, for the matrices in the 
Boltzmann equation we have 



, Mo 

\Mi My' 

It proves convenient to define a 2 x 2 matrix 

_ f Mo -Bo 
[ Bo Mo 

as well as 2 X (27V) and (27V) x 2 matrices TZi,TZi 



Bo Bf 
Bi B 



7^l 



B{ M^ 



7^l = 



Afi -Si 
Bi Ml 



(7.2) 



(7.3) 



(7.4) 



This allows us to reformulate the matrix equation, split- 
ting it into a zero component. 




and a vector component, 




Defining 



g = 



d 



M -B 







7^l 



B M 

we can formally solve for the vector 



= Q 







-7^l 



1 



(7.5) 



1,1 . (7.6) 



(7.7) 



{7- 



Inserting in (j7.5p and solving for the zero mode compo- 
nents leaves us with 







The response functions can finally be calculated from the 
solution for V'ji'i- 
A. Magnetotransport in the absence of disorder 

Significant progress can be made due to a great sim- 
plification which occurs in the above formulae if there is 
no disorder. As we have shown in (|4.29p . in this case the 
vectors Mi and Bi only have one non-vanishing compo- 
nent. Ml = Ml 6*1 and Bi = Biei. The same holds in 
general for the driving terms, = F['ei. One can then 
easily convince oneself that one only needs to know the 
2x2 matrix 



G 



^ ^1 



^ )S{ei ei ) = 



- I 91 92 
'92 91 



(7.10) 



where 



5i(w,B) = ef-(M + BM-iB) ^ei, 
52(w,B) = • M-^B (M + BM-^B) 



(7.11) 



ei, 



and the 2x2 matrix 



Ri 




(7.12) 



both acting in (|| , _L)-space. has been defined in 

App. |B1 Further note that for B — 0, gi coincides with 
the transport coefficient defined in the previous sections. 
For convenience we give the explicit form of i?o (|7.3p us- 
ing results from App. |BJ 



Ro 



2 

N 



-iuj{e + P) 
Bp 



-Bp 

'iuj{£ + P) 



(7.13) 



The above equations immediately yield the solution for 
the and 1 components of -i/jn j^. 



^11 



f) 



— [i?o — RiGR]\ 



G 



pD 







— RiG 



Ri 



(7.14) 



(7.9) 



It is a simple matter to express all response coefficients 
using (|4.22p . For the thermal conductivity, e.g., one finds 
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N 



[(i?.o — RiGRi) ^i?iG] I 



° ^' \-[GRi{Ra-RiGRi)-^] [G + GRi{Rq ~ RiGRi)'^Rig] 




N 



[(i?o — RiGRi) ^ 
[G-Ri (i?o — R\GRi 



[G + Gi?i(i?o - i?iGi?i)-ii?ib] 



.(7.15) 



These expressions are further analyzed below. 



B. D.C. response of pure samples in a magnetic 
field 



1. Small fields - Hydrodynamic regime 



The hydrodynamic regime corresponds to the limit 
where one may approximate 



91 

92 



0- 



(7.16) 



We will see that corrections to higher order in the mag- 
netic field correspond to corrections beyond the hydrody- 
namic analysis, as the latter indeed relies on the smallness 
of the magnetic field. In Eq. (|7.34p below we will give a 
quantitative criterion for the onset of corrections in large 
fields. 

At finite magnetic field and in the absence of disorder, 
the Boltzmann transport theory predicts a vanishing d.c. 
conductivity and thermopower, 



cTxxii^ = 0) = a^xii^ = 0) = 0. 



(7.17) 



This can be understood as a consequence of Lorentz in- 
variance: In a reference frame moving at the constant ve- 
locity vu = E X B/_B^ with respect to the lab frame, the 
observed electric field vanishes, and hence, in that frame 
all currents vanish. Upon transforming back, since v/j is 
perpendicular to E, this implies also a vanishing longitu- 
dinal response to the electric field in the lab frame. The 
transverse d.c. response takes an equally simple form: It 
yields the standard Hall effect, axy = p/B, and the trans- 
verse Peltier effect axy = s/B, which can be interpreted 
as charge and entropy density drifting with the velocity 
V u — . These results are in agreement with the hydrody- 
namic description of Sec. Illli but hold much more gener- 
ally due to Lorentz invariance, even when lu^^'^Tcc 2> 1. 

In contrast to the electrical conductivity and the ther- 
mopower, the longitudinal thermal conductivity of clean, 
interacting Dirac particles remains finite in the d.c. limit. 
From the formula (|7.15p and the matrices (|7.12l7.13p one 
finds after some algebra the result 



The expression simplifies further at particle hole symme- 
try where p = and the matrix G is diagonal (52 = 0). 
This is a consequence of M and B being symmetric and 
antisymmetric with respect to A — > — A, respectively. g2 
vanishes since it is an expectation value of a matrix under 
antisymmetric under A ^ — A, evaluated on the mode 0i 
which has definite symmetry under the same transforma- 
tion. 

Using agifi = 0) = gi[I^+'']X=o = O.lSeeVto^ (for 
1] = 0.01), we find 



Kxx{t^ = 0, = 0) = Kxx{i^ = 0,p. = 0) (7.18) 



1 s'T 



B^ ^q(O) 



= Gmhoq;^ ( —pTT I ^bF, 



\eBvj 



with 



G 



1 



MHD — 



0.156 \ n J 



76.0. 



(7.19) 



Note that the expression (|7.18p is precisely the prediction 
of magnetohydrodynamics, cf. Eq. (|3.1ip . 

The transverse response at particle-hole symmetry is 
found to vanish in a closed circuit. 



,(w = 0,M^O) = 0, 



(7.20) 



in agreement with hydrodynamics as well. However, Kxy 
diverges, since 



Kxyi^^ = 0, ^ = 0) 



'xy 



B 
P 



(7.21) 



Note that, similarly as in the hydrodynamic approach, 
the limits B ^ and a; — )■ do not commute. That 
is, the small field limit of the d.c. response discussed 
above does not correspond to the d.c. limit of the B ^ 
response. Indeed, the latter diverges in general as 1/lu 
due to momentum conservation except at zero doping, 
whereas any finite field B leads to a vanishing d.c. con- 
ductivity in clean systems. 



2. Large fields - Ballistic regime 

The Boltzmann approach allows us to go beyond the 
hydrodynamic regime. In particular, it is interesting to 
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study the large field limit, and how the crossover, as con- 
trolled by the parameter hjo?, takes place. The static 
response in the regime where this parameter is small and 
the inelastic scattering time is the shortest timescale, was 
discussed above. Further, in the subsequent section we 
will show that the complete low frequency response is well 
described by relativistic magnetohydrodynamics, too. 

The regime of large fields can be addressed by deter- 
mining the deviations from equilibrium in a perturba- 
tive expansion in o? jh. Such an approach was taken in 
Ref. llJ where the thermoelectric response of relativistic 
bosons was studied. Those are the low energy quasipar- 
ticles expected at the transition from a commensurate 
Mott insulator to a superfluid, for which the simplest 
model is the Wilson Fisher fixed point occurring in a 0^ 
theory. In an e = 3 — d expansion^^'^^ around d = 3, the 
fixed point value of the quartic coupling scales as e, which 
takes the role of a in our case. In the spirit of an e expan- 
sion, the authors of Ref. [Isl assumed e to be small and 
thus focussed on the 'ballistic' limit of magnetoresponse 
where the parameter is large. Here we analyze the 
large field limit in an analogous way, and work out the 
crossover to the hydrodynamic low field regime. 

We recall that the field should still be weak enough to 
avoid the effects of Landau quantization. This requires 
that the typical cyclotron energy, Sw'^p be smaller than 
fcsT, or equivalently, & <C 1. 

In the regime, ^ & ^ 1 the transverse deviation 
of the distribution functions leads to a transverse cur- 
rent whose deflection by the magnetic field balances the 
current-driving tendency of the applied fields. The lon- 
gitudinal currents are smaller by an extra factor of o? jh. 
Formally, to leading order, the solution of the Boltzmann 
equation (|4.4p is given by 



^11 



(7.22) 



B-^M^'^B-^T. (7.23) 



Let us focus on the thermal response again. Writing 
the driving thermal gradient as 



•^(fc)--/Afc[l-/A\]5T(A,fc)VT.ek 
where 5t(A, k) — k — Xfi we find from (|7.22p 



(7.24) 




400 a 



Figure 5: The thermal conductivity at particle-hole symme- 
try scales B ^ both at small and large B. The 
plot shows the coefficient C in the relation ^^^{fj, — 0; B) = 
C o? /b^k^T as a function of b/a^ ~ u/'J^Tac- It interpo- 
lates between the magnetohydrodynamic regime (|7.19p and 
the large field hmit of Eq. (E28j. 



with a different numerical prefactor: 



K:r2;(t^ = 0, ^ = 0) = Kjja; (tJ 0, ^ = 0) (7.27) 



eBvl 



klT. 



Here, 



(7.28) 



N 



E 



dpk 



2 i2n) 



g(A, fc)X,^,UA,fc I {.g}] = 383.98, 



where g(A, fc) = Afc^, and we have regularized the colli- 
sion term by fixing the RPA screening to rj = 0.01. The 
crossover from the magnetohydrodynamic regime to the 
regime of large fields has been calculated using the full 
matrix formalism of the previous section with 12 basis 
functions, and the result is plotted in Fig.O Notice that 
the crossover extends over rather large fields before the 
asymptotic limit (|7.28p is reached. 



\k 

3±(A,fc) = -^.gT(A,fc)|VT| 



(7.25) 



The longitudinal component follows form the balance be- 
tween the defiection of the longitudinal current and the 
inelastic collisions which degrade the transverse currents 

\k 

g||(A,fc)/,",[l - = - - [Al^V] (A,fc). (7.26) 

We focus on the particle-hole symmetric case where from 
(|7:25l) g±{X,k) \P{-VT)/b. The longitudinal ther- 
mal conductivity k^x turns out to scale parametrically in 
the same way as in the hydrodynamic regime, however 



C. A.C. response in a magnetic field 



Cyclotron resonance 



From the expression (|7.15p it is easy to infer that all 
response functions have a pole in the complex frequency 
plane where the determinant 



Det[i?o - RiGRi] = 



(7.29) 



vanishes. This defines the cyclotron resonance lo* as a 
function of B and p. Its trace in the response is illus- 
trated in Fig. [S] where real and imaginary parts of the 
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Re[axJ,Iml0xx] 




Figure 6: Collective cyclotron resonance evaluated for a small 
magnetic field hjc? = 0.5 and in the quantum critical regime 
^/T = 1 (light doping). We plot the real and imaginary part 
of axx{i^) in units of e^/ft as a function of frequency. The red 
curves were obtained from the full solution of the Boltzmann 
equation, while the blue curves are the magnetohydrodynamic 
prediction. The latter relies on the smallness of which 
is seen to be an excellent approximation for the parameter 
chosen here. The maximum in the real part occurs at the col- 
lective cyclotron frequency Uc = pB/{e + P). The collisions 
introduce an intrinsic damping and lead to a broadening of 
the resonance scaling like at small enough fields. 



longitudinal conductivity are plotted as a function of fre- 
quency, exhibiting a clear resonance. 

The above condition (|7.29p is equivalent to 



e + P 



(7.30) 

2 



At small B the solution is given by 



v(0) 



with 



e + P 



7("'=aQ(/i,0) 



B^ 



e + P 



(7.31) 

0{B'), 

(7.32) 



where aQ{iJ,,uj = 0) was given in IjS.Sp . Note that the ex- 
pression to order 0(3"^) is in precise agreement with the 
predictions from magnetohydrodynamics, cf. Eq. (j3.13|) . 
However, here we have the additional benefit of obtaining 
precise information on the dependence of the broadening 
7 ~ ctq on the chemical potential, as well as on the lead- 
ing corrections in large fields. 

The evolution of the cyclotron resonance with increas- 
ing magnetic field is shown in Fig. [7] for three different 
values of small fields. Note that the peak value of the 
resonance decreases as while its width increases as 

B\ 



Figure 7; Real and imaginary part of a^xii^) (in units of e^/fi 
for fi/T = 1 in magnetic fields b/a^ = 0.75, 1, 1.5 (from left to 
right). The cyclotron resonance ujc is essentially proportional 
to B. The pair of red curves associated to each field value 
corresponds to the full solution of the Boltzmann equation. 
The blue curves are the prediction of magnetohydrodynamics 
which is seen to be an excellent approximation at these small 
field strengths. 



We can also include the effect of weak disorder in the 
determination of the cyclotron resonance. In the most 
general case we would have to look for a zero of 



Det[7^o - niGTZi] = 



(7.33) 



in the complex frequency plane. 

However, similarly as in the calculation of the impu- 
rity scattering rate (|5.15p . one can check that to 
lowest order in A, only the impurity scattering of the 
momentum mode enters via the matrix element M™'', 
and we can restrict ourselves to the analysis of (|7.29p 
with modified Mq. Including the latter into the equa- 
tion for w*, one finds that the pole is simplify shifted by 
«Tj~jp. Such an effect was anticipated by the inclusion 
of a phenomenological momentum relaxation rate in the 
relativistic hydrodynamic o ^ ^ 1 , and is put on a rigorous 
basis here. Of course, the above discussion of disorder 
effects is meaningful only as long as the impurity scat- 
tering rate is the smallest scale in the problem, i.e., for 
A < fe, a^. 

To avoid confusion, we point out that the collective 
resonance frequency ojc differs from w'^p (|3.3p , which we 
had defined as the cyclotron frequency of non-interacting 
particles at thermal energies. The difference is partic- 
ularly marked in the critical regime /i ^ T where tOc 
is significantly smaller than the cyclotron frequency of 
non-interacting thermal particles. It describes the orbit- 
ing motion of the strongly colliding relativistic plasma as 
a whole, where single particles do not have enough time 
to complete an orbit, but constantly collide with others, 
undergoing a collective motion with average frequency 
ujc P proportional to the doped density. 
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B 



2 4 6 8 «^ 

Figure 8: Real part of the cyclotron pole uj* = Uc — i'y, in units 
of ksT /ho? , as a function of magnetic field for fixed chemical 
potential ^/T — 1. The low field part is correctly predicted 
by relativistic magnetohydrodynamics, while at high fields ujc 
tends to exceed a;^^° 



2. Large fields - beyond hydrodynamics 

Inspecting (|7.3ip . an estimate for the field strength 
where higher order corrections in B become important 
can be obtained from either of the conditions 

T'^'gi/i^' ^ 1, or, y")^^ (7.34) 

9i 

One can check that both are equivalent to ^ 1, 

which in turn is equivalent to the physically intuitive con- 
dition 

uj^^Teo - 1, (7.35) 

both in the quantum critical and the Fermi liquid regime. 

In order to characterize the large field regime, we have 
determined the precise location of the cyclotron reso- 
nance by solving numerically for the zero of the deter- 
minant (|7.29p for fi/T = 1, relying again on 12 basis 
functions to calculate gi^2{oJ,B). The result is plotted 
in Figs. IVII C 21 and IVII C 21 We can interpret the in- 
crease of LUc beyond the hydrodynamic expression as an 
effect of oj^^ becoming of the order of t^^: The larger 
the field, the further single particles proceed in their un- 
perturbed, non-interacting cyclotron orbits. Since the 
cyclotron frequency of non-interacting thermal particles 
is typically higher than Wc, one should expect that the 
decrease of the scattering probability leads to an increase 
of the observed cyclotron resonance frequency, as is seen 
in Fig. IVII C 21 Given that the scattering is not very effi- 
cient in this large field regime, it is also natural that the 
hydrodynamic approach overestimates the broadening of 
the resonance due to inelastic collisions, cf. Fig. IVII C 21 

As we pointed out, the magnetohydrodynamic predic- 
tions break down at large fields. This is illustrated in 



y 




2 4 6 8 



Figure 9: Imaginary part of the cyclotron pole uj* = luc ~ i'y, 
in units of ksT/ha^, as a function of magnetic field for fixed 
chemical potential ^/T = 1. The damping at small fields 
is correctly predicted by relativistic magnetohydrodynamics, 
while it is overestimated at large values of b/a^. 



Re[axJ,Im[OxJ 




Figure 10: Real and imaginary part of axx{i^) in units of e^/ft 
in a large magnetic field b/a'^ — 10, at /i/T = 1. The pair of 
red curves correspond to the full solution of the Boltzmann 
equation. The blue curves are the prediction of magneto- 
hydrodynamics which fails completely for this large value of 
b/a^ > 1. 

Fig. [To] for the cyclotron resonance in the conductivity 
which is significantly shifted with respect to the hydro- 
dynamic prediction. 

The result (|7.3ip for the cyclotron pole is generally 
valid at small fields for any relativistic fiuid, as was shown 
by the relativistic hydrodynamic analysis in Refs. flllll^ . 
Quite remarkably, this cyclotron resonance also emerged 
from the solution of an exactly soluble, but strongly 
coupled supersymmetric conformal field theory via the 
AdS-CFT mappingi^. Furthermore, the deviations from 
the hydrodynamic prediction at high fields could be 
found numerically in that case as well, and the trends 
of ll!c{B), j{B) were found to be very similar to those in 
Figs. IVII C 2|VII C 21 This is very interesting since in the 
present work we are limited to weak coupling a <C 1, for 
the Boltzmann approach to be quantitatively accurate. 

We note that the cyclotron resonance should be readily 
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Figure 11: Sharp cyclotron resonance in the conductivity (in 
units of e^/h) for fi/T = 5 and a moderate magnetic field 
b/a^ — 1. Since in the Fermi liquid regime the presence of 
antiparticles (holes) can be neglected, and the relevant part 
of the single particle dispersion is hard to distinguish from a 
parabolic dispersion, Kohn's theorem applies approximately, 
the cyclotron resonance becoming sharper and sharper at 
asymptotically large /i/T. 



observable in graphene at T of the order of room tem- 
perature and in moderate magnetic fields corresponding 
to fractions of a Tesla, as discussed in Ref. 



D. Recovery of Kohn's theorem 

The above collective cyclotron effects are most pro- 
nounced in the relativistic, quantum critical regime < 
T. As one leaves the latter, the cyclotron resonance at a 
given magnetic field B becomes sharper and sharper as 
characterized by the ratio which equals the number 
of applied flux quanta per doped carrier, multiplied by 
CQ^fi) (as measured in units of e^/h). As we have seen 
above, (Jq{ij) decays as {T/^)"^ at large fi. The result- 
ing sharpening of the resonance is clearly seen in Fig. 1111 
which shows axxi^) evaluated for fi/T = 5. 

It is interesting to note that as fi/T increases, the res- 
onance approaches the value 



,(0) = 



eB 



pB eB 
: + P P'/v'^p hkp/vp' 



(7.36) 



which one recognizes as the semiclassical cyclotron fre- 
quency expected for a circular Fermi surface at wavevec- 
tor kp and Fermi velocity vp. Both observations indicate 
that one recovers the familiar Fermi liquid characteristics 
at large doping. 

The intrinsic broadening of the cyclotron resonance 
due to collisions is an interesting effect pertaining mostly 
to the quantum critical regime < T. As noted above, 
in the Fermi liquid regime the width 7 of the cyclotron 
resonance tends to zero. This can be understood as re- 
flecting the approach of a regime where Kohn's theorem 
should apply asymptotically: One species of quasiparti- 



cles is entirely frozen out in this regime, and the devia- 
tion of the linear band structure from a parabolic disper- 
sion becomes increasingly negligible. These are the two 
conditions under which Kohn's theorem is valid: The 
latter asserts that in a system with only one parabolic 
band there is a single sharp resonance peak at a well de- 
fined cyclotron frequency, irrespective of the presence of 
electron-electron interactions. 



E. 



Range of vahdity of relativistic 
magnetohydrodynamics 



Using Eqs. (|7.14l7.15p one can obtain explicit expres- 
sions for the frequency dependent response in pure sys- 
tems in a magnetic field. Note the remarkable fact that 
the three pairs of longitudinal and transverse response 
coefficients are strongly constrained: Apart from ther- 
modynamic data such as {s + P), s,T, p and fj, the for- 
mulae contain only two independent matrix elements 
(7i^2(w, B, /i). To leading order in magnetic field one can 
even neglect the dependence on 52, and all the response 
functions arc interrelated, with one single parameter left, 
corresponding directly to erg in the hydrodynamic formu- 
lation. 

The full expressions for the response functions can eas- 
ily be worked out analytically, but the expressions are rel- 
atively lengthy and will not be displayed here. However, 
it is interesting to use these exact results to determine 
the extent to which the magnetohydrodynamic formu- 
lae given in Sec. IIIII are valid. As we already know, at 
large fields corrections set in, and similar corrections are 
to be expected at higher frequencies of the order of the 
inelastic scattering rate. We are now in the position to 
characterize the corrections precisely. A detailed anal- 
ysis of the response shows that the following statement 
holds: For small frequencies {yjjcP' ^ 1) and small fields 
(6/a^ ^ 1), and a fixed ratio the exact a.c. conduc- 
tivity satisfies the asymptotic equality 



,(w,B) 



.MHD 



{uj,B)+0{hlo?,i^lo?). (7.37) 



An analogous relation holds for all other response func- 
tions. Hereby the expressions (|7.32|1 . with oq from (|5.5p . 
have to be used in the magnetohydrodynamic response 
functions given in Sec. IIIII 

It is interesting that a very similar result was obtained 
from the exact solution of the strongly coupled conformal 
field theory studied in Ref. [IB, showing that the validity 
of this statement is not restricted to weak coupling. 



VIII. CONCLUSION 

We have presented a Boltzmann approach to describe 
transport in liquids of interacting Dirac fermions with 
and without magnetic fields. We have established that 
as long as the inelastic scattering rate is the largest 
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scattering rate in the problem, the relativistic hydrody- 
namic formalism captures the frequency-dependent re- 
sponse very well. Further, we have obtained an exact 
expression for the single transport coefficient uq that is 
left undetermined by hydrodynamics, and showed that it 
decays as a power law as one leaves the quantum critical 
relativistic regime. At large doping the electron system 
was shown to recover all the signatures of a Fermi liquid, 
such as Mott's law and the Wiedemann Franz relation. 
At the same time the collective cyclotron resonance, a 
remarkable feature of quantum criticality, turns gradu- 
ally into a sharp resonance centered at the semiclassical 
cyclotron frequency as one dopes the system further. 

Finally an analysis of the large field behavior yielded 
similar qualitative deviations of the cyclotron pole from 
the corresponding hydrodynamic prediction as was found 
in the exact solution of a strongly coupled conformal field 
theory. 
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role of the speed of light is taken by the Fermi velocity 
vp. The velocity field (in units oi vp) is determined 
in such a way that there is no energy flow in the local rest 
frame where u'^ = (1, 0). Notice that due to the presence 
of heat flows this does not coincide in general with the 
velocity defined by the average charge current. 

The additional terms in (|AllA2p are dissipative con- 
tributions: the Reynolds tensor r^'' accounts for viscous 
forces, which turn out to be irrelevant for small wavevec- 
tor response. The vector v^^ is proportional to the heat 
current. To obtain a closed set of equations, one has to 
express the heat current in terms of local quantities. In 
a relativistic system, its form is strongly constrained by 
covariance and the requirement that the entropy of the 
liquid always increases. The divergence of the entropy 
current follows from the equations of motion as: 



SU V 

T 



(f ) - lfF„ 



T 



(A5) 



where we have used the thermodynamic identity e + P — 
fj-p = sT for the entropy density, s. The requirement 
that the left hand side be positive, and the assumption 
that the heat current should be linear in the gradients 
of T and as well as in the electromagnetic fields (i.e., 
the gradients of the scalar and vector potential), imposes 
that 



Appendix A: RELATIVISTIC HYDRODYNAMICS 

In this appendix we discuss the hydrodynamic and con- 
stitutive equations which are used to obtain the hydro- 
dynamic response functions given in the main part of the 
text. 

In covariant notation, the conservation laws for a rel- 
ativistic fluid read 



9„J" = 0, 



(Al) 

(A2) 



where the energy-momentum tensor and current vector 
of the fluid are given by 



(A3) 
(A4) 



where e is the energy density, P is the pressure, p the 
charge density, g^'^ = diag(— 1,1,1) the Lorentz metric 
and F^'' is the electromagnetic field tensor. Note that the 



^" =-<TQ[g''^ +u''u^)[Tdx{^)-F),^u<]. (A6) 

A similar relation holds for r"'''. This leaves us with a 
single undetermined transport coefficient erg > with 
units of a conductivity. Note that in the relativistic case 
the heat current is not only proportional to the thermal 
gradient, but also to the acceleration of the fiuid element 
(second term). Note that the assumption that v"" is linear 
in F^'' restricts the above argument to small fields B. 



Appendix B: MATRIX ELEMENTS OF TERMS IN 
THE BOLTZMANN EQUATION 

1. Matrix elements in a general basis 

In this appendix we give explicit expressions for the 
matrix elements appearing in the Boltzmann equation 
after a projection onto a specific basis of modes. The 
collision integral describing electron-electron scattering 
as well as impurity scattering is given by 
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Icon[A,A,t|{/}] = 2vra^mLA^| (Bl) 

5{k - fci - |k + q| + |ki - q|)i?i(k, ki, q){/A(k, i)/_A(ki, i)[l - /A(k + q, t)] 

X [1 - /_;,(ki - q, <)] - [1 - A(k, t)][l - /_A(ki, i)]/A(k + q, i)/-A(ki - q, <)} 
<5(fc + fci - |k + q| - |ki -q|)i?2(k,ki,q){/A(k,i)/A(ki,i)[l -/A(k + q,t)] 

x[l - /A(ki - q,t)] - [1 - /A(k,t)][l - /A(ki,t)]/A(k + q,t)/A(ki - q,t)}| 
+27r 1 ^<5(fc-fci)|C/AAp[/A(k,t)(l-/A(ki,i))-(l-/A(k,i))/A(ki,t)' 



where 



i?i(k,ki,q) = — (|r+__+(k,ki,q)-r+_+_(k,ki,-k-q + ki)|' 

+ (7V-l)|T+__+(k,ki,q)f + (7V-l)|r+_+_(k,ki,-k-q + ki)|') , 



4 /1| i^i 
i?2(k,ki,q) = ^ - r++++(k,ki,q) - T++++(k, ki, ki - k - q) 

+ (iV-l)|T++++(k,ki,q)|') , (B2) 



and the disorder potential ?7aa was introduced in Eq. (I2.13p . The inelastic scattering rate is proportional to a^, and 
the matrix entries of the electron-electron scattering matrix A^*^*^ read explicitly 



^ inn 



2-"' E / (01^ (l^-^^^ + fci - |k + q| - |ki - q|)0„(A, k) (B3) 
|i?i(k,q - ki,q)(l - fl,){l - /°A|k,-q|)/-Afe,/A%+q| 

xek • [ek(/>«(A, fc) + eki0ri(-A, fci) - e]^+q^ct>n{\ |k + q|) - eki-q0„(-A, |ki - q|)] 

+i?2(k,ki,q)(l - /°,)(1 - /A\)/AV+q|/AVl-q| 

XGk • [ek0„(A, k) + eki(?!)„(A, fci) - ek+q(?!)„(A, |k + q|) - eki-q0„(A, |ki - q|)]| , 

I 



where The matrix for the time derivative reads 

,U ■ ^ ; ■ A^imp -^™„=*^E / 7TT2/Afe(l-/Afe)<^™(A,fc)0«(A,fc) 

Ihe entries m the impurity scattering matrix y\4 ^ a^ l^"') 

are given by 



(B7) 

'21, ^2; 



(27r) (27r) |k ki| Finally, the deflection of currents by the magnetic field 

x/Afc(l - /AfcJ'^m(A, k) [0„(A, fc) - Gk • eki0„(A, fci)] , is described by a matrix with entries 
where the strength of disorder is measured by the dimen- 

sionless parameter ^^^^^^ ^ f <fk^,^^^ /a",)^0„.(A, fc)0„(A, fc) , 

"^^"Mfc^j (B8) 
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where the dimensionless parameter measuring the field 
strength is given by 



6= 



(B9) 



The projection of the driving terms due to an electric 
field or a temperature gradient onto the basis functions 
yields the vectors 



J- — 



A 

-E 



.fLi^-fxM\k), (BIO) 

:(fc-AM)/O,(l-/O,)0™(A,fc). 



2. Two mode approximation 

Here, we evaluate the matrix elements with respect to 
the two main modes, (j>o,(j>i, of the specific basis (|4.25t - 
I4.27p . In the sector spanned by these two modes, the 
electron-electron collision operator takes the form 



M 



Cb 




Mf^(^) 



(Bll) 



where 



Mfn.)^«^ETy(2^(2^(2^{ 
5{k - fci - |k + q| + |ki - q|)i?i(k, ki, q) 
^(1 ^ /"*;)(! ^ /-Afci)/A|k+q|/-A|ki-q| 
X "k ki (k + q) ^ (ki-q) l^ 

[k ki |k + q| |ki-q|. 
+Sik + fci - |k + q| - |ki - q|)i?2(k, ki, q) 
X(l — /Afe)(l ~ /Afei)/A|k+q|/A|ki-q| 

\k ^ ki (k + q) (ki 

k ki 



|k + q| 



|kl 



(B12) 



The matrix elements of all other operators can be ex- 
pressed with the help of the functions 



N 

y 



dk 



i6s^+ + Xds^-)k''l 



n 

\k 



2n ^ J 2tt 

where a partial integration was used to obtain the second 



line. One easily verifies the explicit relations 



7-(0) _ -• - 



j(0) ^ 



7-(2) _ 

7-(2) _ 



N 1 
2" 2^' 
N tanh(/i/2) 

Y 2^ ' 
N 1 

--ln[2cosh(^/2)], 

Z TT 

N [I 
Y2^' 

+ 



(3) 



- p- ^ p 

2 ' 



(B13) 
(B14) 
(B15) 

(B16) 

(B17) 
(B18) 

(B19) 



where p^ are the number densities of particles and holes, 
respectively, and p is the charge density in units of e. The 
last relation follows since in a two-dimensional relativistic 
liquid P — e/2. This is a consequence of the energy 
momentum tensor being traceless, a relation that can 
easily be checked explicitly for free Dirac fermions. 

In order to analyze the Fermi-liquid regime it is conve- 
nient to have the asymptotic form of /j^ for large p at 
hand. From a standard Sommerfeld expansion one finds 



/f^(M) 



N 1 

'2 2^ 

N p 

'2 2^ 



N 
N 



2tt 
27r 



2^' 



(B20) 



up to corrections of order 0(exp(— /^)). 

In the two mode approximation, the scattering from 
Coulomb impurities is described by 



/f (M) /i^^(M) 



(B21) 



and the matrix for the time derivative takes the form 



N 



e + P p 



(B22) 



Finally, the matrix describing the deflection by the mag- 
netic field has the two mode representation 



B^b— 



N\ll'\p) I'^^ip) 



0)l 



(B23) 



Due to the choice (|4.27p the driving fields only have 
components along c/jq and given by 



P 

N\l^^\p) 



(B24) 
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and, 



2 / e + P-np 
N f ¥ln[2cosh(^/2)] 



(B25) 



Appendix C: INELASTIC SCATTERING RATE IN 
THE FERMI LIQUID REGIME 

In this appendix we analyze the behavior of the ma- 
trix element in the limit of large chemical potential 
H ^ T, in order to estimate the quantity g^^ which de- 
termines the inelastic scattering rate. 

The phase space for the outer k-integral in (jB12p scales 
like p,T. Further, given a momentum transfer q, the inte- 
gral over ki over the energy conserving 5-function scales 
like Tfi/\q\, where the factor l/|q| is due to the derivative 
of the argument of the (S-function. The factors Ri,2 scale 
like a^/lqp if the Coulomb interactions are not screened, 
and thus, for small q the integral over the momentum 
transfer behaves like / (Pqq~^ . The infrared divergence 
is cut off by the square of the difference of distribution 
functions g\{k) which provides an extra factor of at 
small q. The cut off scale is set by the typical range over 
which g{k -\- q) varies, which is 9 ~ T. Consequently the 
integral is dominated by g ~ T which contributes a phase 
space factor . In the final estimate of gi^{^) we need 
to take into account that the relevant mode for electrical 
conductivity, (pi , has a strong overlap with the zero mode 
(/)o , when k is restricted to the thermally relevant vicinity 
of the Fermi level ji. Since is a zero mode of the in- 
tegral, only the part of (pi orthogonal to (^0 contributes, 
and this provides an extra factor of [T / . Putting all 
these factors together and multiplying by a normalizing 



factor 1/T'^, we find that unscreened electron-electron in- 
teractions are dominated by small momentum transfer of 
order g ~ T, leading to a scattering parameter 



9i ^ip) - ^ 



q |q|2 



= (CI) 



which tends to a constant. 

The inelastic scattering rate can be estimated from a 
similar expression, where the first k-integral and the nor- 
malization factor are dropped. Looking at the relaxation 
of modes g{k) that have a variation of order 0(1) over an 
interval of order 0{T) around k — fj,, there is no suppres- 
sion factor from the square of differences in g\{k)^ and 
we obtain the estimate 



Tpq^ 

1 |q|' 



(C2) 



This is indeed consistent with the expression (|5.6p . upon 
using the estimate (|5.7p . 

If we include screening of the interactions, the scatter- 
ing is eventually dominated by q ~ /i, while the phase 
space integral of the ki -integral contributes qT instead 
of q^. This leads to 



q |q|2 ^2 



9~M 



= 0(a^^),(C3) 



and the analogous estimate for Tcc yields the familiar 
Fermi liquid behavior 



a 



TjiqT_ 



(C4) 
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